Thermal conductivity in the vortex state of YBa2Cu307 and Sr2Ru04 
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The thermal conductivity k is calculated for the Abrikosov vortex lattice with d-wave pairing. 
The Andreev scattering rate arises from the self energy term in Gorkov's integral equation for the 
I corresponding Green's function. The impurity scattering rate is calculated self-consistently from 

this Green's function in the t-matrix approximation. Neither point-like Impurity scattering with a 
' phase shift in the almost unitary limit nor in the Born limit can consistently explain recent zero 

field and finite field measurements of k/T at low T in ultraclean YBCO. Our theory also applies to 
the spin triplet pairing states with vertical or horizontal line nodes in Sr2Ru04 . 
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5^ , Recently the thermal conductivity k of ultraclean YBa2Cu307 wa.s measured at very low temperatures in magnetic 
rS i' fields up to 13 Tesla.^ In zero field, k/T rises from its universal limit koq/T very fast as a function of temperature 
^ ' which, according to the quasiclassical calculation for a d-wave superconductor, ^ indicates either an extremely small 
normal state impurity scattering rate F = 1/2t in the unitary limit (phase shift p = 7r/2), or an unrealistically large F 

■ for the Born approximation (phase shift p — ). As a function of field, n initially increases very rapidly reaching nearly 

a twice its universal limit kqo and remains almost unchanged up to 13 Tesla. The theory of Vekhter and Houghton,'^ 
I ' which includes scattering of quasiparticles at the vortices via Andreev scattering, yields a rise proportional to -v/ff in 
, agreement with the measurements for a less pure sample of YBCO. ^ Recently the data of Ref. 1 have been accounted 
' for by assuming a phase shift p for impurity scattering slightly less than the unitary limit, 7r/2 , and adding to the 
O I extremely small impurity scattering rate F a phenomenological quasiparticle scattering due to Andreev scattering at 
I the superfliud flow of the vorticies.^ 

In the present paper we use the theory of Ref. 3 and the equivalent theory of Ref. 6 to solve the problem of the 

■ occurence of a plateau-like feature in the field dependence of k in YBCO. These theories are based on the theory 
' of Brandt, Pesch, and Tewordt (BPT-theory) ^ where the spatial average of the Gorkov Green's function G for the 

Abrikosov vortex lattice was calculated. The expression for the thermal conductivity which has been derived from this 
• Green's function^ can explain the measured field dependence of k for different field directions for superconducting 
I ' spin-triplet pairing states with vertical or horizontal nodes in Sr2Ru04 . A simplified version of the BPT-theory which 
is based on the quasiclassical Green's function g , referred to here as the P-approximation, yields the expression 
of Ref. 3 for the thermal conductivity the vortex state with the field perpendicular to the basel plane for 

a d-wave superconductor. We find that the expressions for the thermal conductivity derived from the BPT-Gorkov 
or from the P-quasiclassical Green's functions (Refs. 6 and 3) yield to a very good approximation the same results 
^ for all reduced fields h = H/Hc2 , impurity scattering rates S = F/Ao , and phase shifts p. Here Aq is the amplitude 
of the d-wave order parameter. 

To save space we omit here the full expression for k and concentrate on the discussion of the Andreev scattering 
which is derived from first principles. This contribution to the quasiparticle scattering is contained in the expression 
Q ^ for Im^o in the denominator of the w-integral for n which also includes the well-known factor sech^(w/2T) Here 
Q . is the position of the pole of the BPT-Green's function G as a function of the normal state energy ^ measured from 
the Fermi energy. The equation for the zero of the denominator of G yields Im^o = + lA , where 7i is the total 
impurity scattering rate and "/a is the imaginary part of the quasiparticle energy —Tja at • The kernel of Gorkov's 
integral equation for G in the spatial representation yields 

u • 

I]A(ri,r2;L^) = -A(ri)A*(r2)G°(ri-r2;-u;), (1) 

where A(r) is Abrikosov's vortex lattice order parameter. From this expression it is clear that 7^ — — ImEyi is the 
scattering rate for converting a quasiparticle at A*(r2) by Andreev reflection into a quasihole and then back into a 
quasiparticle at A(ri) . The Fourier transform of — 1^2',^) yields^ 

Sa(p, iu) = -iV^A'^{A/vi_)w[iuj + ij, + Cp)A/wx] • (2) 

Here A^ is the spatial average of |A(r)p , A — {2eH)^^^^ is of the order of the vortex lattice constant, v±{p) is the 
Fermi velocity component perpendicular to the applied field H , and w is Dawson's integral. In the limit uj ^ 
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corresponding to the T ^ limit of k/T, one obtains explicit expressions for Im^o = 7i + 7a(Co) and In the 

limit v± ^ , i.e., v{p) \\ H , Eq.(2) tends to the self energy of the BCS Green's function. In the limit A ^ oo , i.e., 
if — > , the expressions for k in the BPT-approximation^ and in the equivalent P-approximation^ tend correctly to 
the expression for k which was first derived in Ref. 10. 

In the following we present results for a d-wave pairing state in a field perpendicular to the basel plane where the 
Abrikosov vortex lattice state is multiplied by cos(2(/)) . The impurity scattering rate 7i is calculated self-consistently 
in the t-matrix approximation for the self energy Ej : 

7i = ReSj ; Sj = T g{u}) [cos^ P + 9^ sin^ p] ^ ; (3) 



g{uj, (t>)=[l- iV^{2AA/vf cos" {2(l))w' {z)] ^'"^ ; (4) 

2 = 2(a; + iEi)A/w; k = (2eR)-^l'^ . (5) 

Here g is the quasi-classical Green's function in the P-approximation,^ and g is the angular average of g . The field 
dependence of A is approximately given by A = AoVl - h , where h = BjUc^ ■ 

In Fig. 1(a) we show the ratio of k^x = k to the normal state conductivity «;„ vs ft- = H/ Hc2 for w = corresponding 
to the limit T ^ , and the density of states N/Nq = Re ^(0) , for reduced normal state scattering rates 5 = T/Aq = 
0.1, 0.01, 0.001, and 0.0001 in the Born approximation (phase shift p = 0). The occurence of a plateau in the low 
field range for all scattering rates is demonstrated more clearly in Fig. 1(b). In Fig. 2(a) wc show our corresponding 
results for and N/Nq for a phase shift p = 0. 49577 which is very close to the unitary limit p = 7r/2 and has 

been used in Ref. 5 to fit the data of Ref. 1. Fig. 2(b) shows more clearly that a plateau in the low field range is 
reached only approximately in the limit of very small scattering rates, here 5 = 0.001 and 0,0001. The latter value is 
of the order of magnitude of the value that has been used in Ref. 5 to fit the data of Ref. 1. We find that the field 
dependence of is very nearly the same as the field dependence of the ratio of the scattering rates in the normal 
and vortex states, i.e., the angular average of V /{^i + 7/i) • This means, according to Figs. 1(a) and 2(a), that the 
Andreev scattering rate is much larger than the impurity scattering F up to fields just below ffc2 , and that the 
ratio 7A/r at fixed /; increases with decreasing scattering rate 6 . 

We now attempt to explain the second feature of the experiments of Ref. 1, i.e., that k/T becomes almost temper- 
ature independent up to about 0.6 K for constant fields above 1 Tesla up to 13 Tesla. In Fig. 1(c) we show our results 
for K,{uj)/Kn vs J7 = w/Aq for 5 = 0.1 and 0.05 and several fields h — 0.1, 0.06, and 0.02 in the Born approximation 
(p = 0). Here, K(a;)/K„ is the factor multiplying (w/T)^ sech^(a;/2T) in the normalized integral over d{uj/T) for 
n/Kn ■ Since (w/T)^ sech^(a;/2T) has a maximum at lo/T = 2.4 , one has Ct = lo/Aq ~ 2.4(T/Ao) , and thus k{u))/k„ 
yields approximately the dependence of k/k„ cx k/T as a function of ~ 2.4(T/Ao) . We have divided k((x')/k„ by 
6 because ~ (7r/2)Koo where kqo is the universal conductivity limit. One recognizes from Fig. 1(c) that K(a;)/K„(5 
is almost constant in the range from Q = to 0.1 which means that k/T is almost constant up to about T/Tc 0.1 , 
and that the curves for constant fields h = 0.1, 0.06, and 0.02 lie close together for each of the two scattering rates 

5 = 0.1 and 0.05. In Fig. 2(c) we show the corresponding results for the phase shift p ~ 0.4957r and scattering rate 

6 = 0.0001 , and in Fig. 3(a) the results for scattering rate 6 = 0.001 and the same phase shift and fields. Again 
these results indicate that k/T is nearly constant as a function of T at least for the higher fields h = 0.1 and 0.06. 
However, these curves for the almost unitary limit are not as close to each other as those for the Born approximation 
in Fig. 1(c). This is in disagreement with the data for the temperature dependence of k/T for fixed fields between 0.8 
and 13 Tesla which lie on about the same level. ^ 

We have also calculated K{u))/Kn6 in the limit of zero field. The expression for k in Ref. 6 yields in the limit H ^ , 
or A ^ oo , the following result: 

- c«»(« 1 ^ i (l + SM^l : (6) 
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n = li + iEi/Ao 
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- cos2(20) 
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n = w/Ao 



(7) 



The impurity self energy is calculated self-consistcntly with the help of Eq.(3). Eq.(6) agrees with the general 
strong-coupling result of Ref. 10. In Fig. 2(c) we have plotted our results for K{iv)/KnS for phase shift p = 0.4957r and 
scattering rate 5 = 0.0001 , and in Fig. 3(a) we show the results for the same phase shift and scattering rates 5 =0.001, 
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0.01, and 0.1. One sees from these figures that for H = the ratio k{u))/ KnS — {2 / tt) k{uj) / kqq tends for O ^ , or 
T , correctly to the vahie (2 /it) , and that this function rises approximately proportional to with a slope that 
decreases for increasing 6 . The zero field and constant field curves in Figs. 2(c) and 3(a) for 6 =0.0001 and 0.001 are 
qualitatively similar to the data for the T dependence of k/T in zero and constant fields.^ However, the values of 
k/T in units of kqq/T shown in Fig. 2(c) for the constant fields h =0.1, 0.06, and 0.02 are an order of magnitude too 
large in comparison to the experimental values. The corresponding values for 5 — 0.001 shown in Fig. 3(a) are much 
smaller. These values can be reduced further to reach the experimental value of about 2 with a reasonable increase in 
the factor a multiplying the rjuantity AA/w [see Eqs. (2) and (4)] for a d-wave pairing state. ^ Here a <x V2/V where 
V2 is the slope of the d-wave gap at the node. The dashed curves in Fig. 3(a) correspond to the value a = 5 . These 
curves are reduced to those in Fig. 3(b) by using the value a = 31 . The upper dashed curves in Fig. 2(c) for 6 = 0.0001 
also correspond to a = 5 and they are reduced to the lower dashes curves by taking the same value of a = 31 as 
for S = 0.001 in Fig. 3(b). A still much larger value of a is needed to reach the experimental value k/koo ^ 2 . That 
corresponds to an unrealistically large value of V2/V , We thus conclude that the measured T dependence of k/T for 
ultraclean YBCO in zero field, and finite fields up to 13 Tesla, can be qualitatively explained by our theory with the 
parameter values p = 0.4957r and 5 = 0.001 and 0.0001. However, the corresponding curves for the field dependence 
of k/ku in the T limit [see Fig. 2(b)] do not rise as abruptly to the plateau-like value for increasing field as the 
experimental curve. Contrary to the nearly unitary limit, in the Born approximation the field dependence of k/«;„ 
shown in Fig. 1(b), and the temperature dependence of k/T [see Fig. 1(c)], are, for example for 5 = 0.1 , in much better 
agreement with the experiments of Ref. 1. However, the T dependence of k/T for zero field is in total disagreement 
with the data: it rises steeply for very small O and then tends to a constant value of about 4koo /T for fl up to about 
0.1 [see Fig. 1(c)]. 

We have also investigated intermediate values of the impurity scattering phase shift, for example, p = OAn which has 

been used to fit the data for the microwave conductivity in YBCO films. ^-'^ We find that the measured field dependence 
of the thermal conductivity^ is better described with p = QAtt and 6 — 0.001 or 0.0001 than with p — 0.4957r and 
and the same 5's. However, the 51 dependence of k/kqo for zero field is quite far from a 17^ dependence: it rises 
steeply to a maximum and then decreases slowly for increasing Q . Thus phase shifts intermediate between the Born 
approximation and the unitary limit do not seem useful for explaining the experiments of Ref. 1. 

The aim of the theory for k in Ref. 6 was to explain the measurements of k in Sr2Ru04 for fields perpendicular to 
the ab plane and for rotating in-plane fields. Recently, k was measured in Sr2Ru04 in zero field down to very low 
temperatures.^^ This low-temperature behaivior of k demonstrates the universal character of the heat transport due 
to a gap with nodes and suggests strong impurity scattering with a phase shift close to tt/2 . Since our present theory 
of K for a d-wave superconductor also applies to a spin triplet f-wave pairing state with vertical or horizontal line 
nodes,^ we conclude that our results for phase shifts p = 0.4957r ~ 7r/2 also apply to Sr2Ru04,. The T dependence 
of k/T for 5 = 0.01 and 0.1 shown in Fig. 3(a) for zero field is similar to the data in Ref. 13. The curves for vs 
H/Hc2 for S = 0.01 ansd 0.1 in Fig. 2(a) lie far below the curves for 6 = 0.5 and 0.2 obtained in the Born approximation 
in Ref. 6. 

In summary, we have calculated the thermal conductivity for the Abrikosov vortex lattice state with d-wave pairing 
which automatically includes Andreev scattering due to the self energy in Gorkov's integral equation for the Green's 
function. The most important result of the present paper is that the BPT-approximation^^^ based on the Gorkov 
equations yields very nearly the same results for k as the P-approximation^'^ based on the equations for the quasiclas- 
sical Green's functions. These results for k correspond to previous results for the density of states N/Nq for d-wave 
pairing in the vortex state. In Ref. 14 it was shown that the BPT- and P-approximations yield nearly the same 
results as the solutions of the quasiclassical equations for applied fields between Hc2 and Hd while the Doppler shift 
approximation gives rise to substantial errors. In agreement with Ref. 3 we find that, for impurity scattering in the 
Born limit, k exhibits a steep initial increase as a function of field and then becomes almost constant (see Fig.l) while, 
in the almost unitary limit, k rises like a/H (see Fig. 2). The latter result disagrees with the results of Ref. 5 where, 
in the almost unitary limit, k rises initially very rapidly and then becomes almost constant with a value in units of 
the universal conductivity koo close to the experimantal value. The reasons for the discrepency arc apparently the 
following. First, the rcnormalization of the impurity scattering self energy in the t-matrix approximation is carried 
out in Ref. 5 by employing the zero-field density of states while, in Ref. 3 and here [see Eqs. (3) and (4)], the density 
of states is calculated self-consistently for finite field. Second, the Andreev scattering rate 7^ is approximated in 
Ref. 5 by bEn where Eh ^ v/A [see Eq.(5)] is the magnetic field energy and 6 is a parameter which is fitted to 
the data of Ref. 1. This expression for the Andreev scattering rate is taken in analogy to the expression derived 
by the method of the Doppler shift followed by an averaging over a single vortex. However, the expression for 
derived in Ref. 15 contains an additional exponential factor which changes the field dependence considerably. In our 
approach "/a = — ImE/i where T,a is the self energy in the Gorkov equation which has the form for Andreev scattering 
by the spatial variation of the complex order parameter A(r) of the total Abrikosov vortex lattice [see Eq.(l)]. It is 
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interesting that the exponential function arising from the imaginary part of Eq.(2) has, in the absence of impurity 
scattering 7^ , a similar form as that of Ref. 15 considered as a function of frequency, field, and angle 9 — Z(p,i?) . 
Our "fA in the absence of impurity scattering is shown as a function of 9 for several different fields and frequencies in 
Fig. 5 of Ref. 16. It should be noted that here 7a depends also on 7^ [see Eq.(2)] where both quantities are calculated 
self-consistently together with Eqs.(3) - (5). Although the theory of Ref. 7 was originally derived for applied fields 
H near Hc2 , it has been shown to work well over the entire range of linear magnetization,^^ and it has been tested 
down to Hci for a d-wave superconductor by comparison with the solutions of the quasiclassical equations.^'' 

Let us now briefly discuss why our theory in its present form fails to explain the experiments of Ref. 1 consistently for 
zero and finite fields. The measured rapid growth of the zero field k/T with T in ultraclean YBCO can, in agreement 
with Ref. 5, indeed be explained only by assuming a very small impurity scattering rate (5 = F/Aq ~ 10""^ to 10~^ 
and a phase shift close to the unitary limit 7r/2 for isotropic scatterers. The observed sudden onset of a "plateau" 
in k/T as a function of field H requires, for this phase shift limit, the assumption of a very small S < 10"'^ because 
this makes the region of V^-behavior small and k beyond this region nearly constant as a function of H . The third 
feature of the experiments, i.e., that k/T is independent of T , is also satisfied. However, the saturation values of k/T 
in units of kqq/T are too large and lie too far apart for different fields in comparison to the experimental values. The 
Born approximation (phase shift zero) with (5 ^ 0.1 to 0.01 yields a much better description of the measured field 
dependence of k , i.e., a steep initial increase followed by a plateau, and the saturation values of k/T in units of kqq/T 
for different fields are in fair agreement with the data. However, the calculated T dependence of k/T in zero field is 
quite different from the observed temperature dependence if one assumes isotropic impurity scattering in the Born 
limit. This is also the case for intermediate phase shifts between p = and p = ■jt/2 , for example for p = OAtt which 
has been used to fit microwave conductivity data.^^ 

In conclusion, we have shown that the measured field dependence of k/T in ultraclean YBCO,^ in particular the 
plateau, can be better described by the Born approximation than by the almost unitary phase shift limit for impurity 
scattering. The use of the almost unitary limit is however necessary ^ in order to explain the observed rapid increase 
oi k/T with T in zero field with the model of point-like impurity scattering with a single phase shift. It is possible 
that the deficiencies of the present theory may be due to this over-simplified model. Similar difficulties in explaining 
microwave conductivity measurements in very clean YBCO have led to consideration of impurity potentials with finite 
range. This gives rise to considerable changes in the density of states at low frequencies due to renormalization of 
the d-wave gap.^^ It seems possible that such effects could lead to a better description of the T dependence of k/T in 
zero field in the Born approximation. On the other hand we believe that the theory of Andreev scattering contained 
in the theory of k in Ref. 3 need not be altered because it agrees with our theory where the Andreev scattering is 
evident from the form of the self energy in the Gorkov equations. 

Our results for the unitary limit and scattering rates 5 0.01 and 0.1 also apply to spin triplet states with vertical 
or horizontal line nodes in Sr2Ru04 . Indeed, the zero field measurments at very low T strongly suggest a phase shift 
close to 7r/2 }^ 

We thank K. Scharnberg for helpful discussions. 
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Fig. 1a 
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FIG. 1. la) Thermal conductivity ratio, k/ku , vs ft = H/Hc2 at T = in the Born approximation for the reduced scattering 
rates S = F/Aq = 0.1, 0.01, 0.001, and 0.0001 (sohd curves, from top to bottom) and the density of states at the Fermi energy, 
N/No vs h (dashed curves from top to bottom). 
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Fig. 1b 
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Fig. 1c 
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FIG. 3. Ic) The factor in the integrand of the a;-integral for k, k{lij)/k„S , vs f2 = uj/Aq . Upper curves for 5 = 0.05 in the 
Born approximation for h = 0.1, 0.06, and 0.02 and lower curves for 6 = 0.1 and the same h (from top to bottom). Solid curve 
for /i = and 6 = 0.1. 
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Fig. 2a 




Fig. 2b 




Fig. 2c 




Fig. 3a 




FIG. 7. 3a) Factor in the w-integral for n , K{ijf)/ n„5 , vs = oi/Ao , for phase shift p = 0.4957r . Sohd curves for ft = and 
S = 0.001, 0.01, and 0.1 (from left to right). Dashed curves for 5 = 0.001 , a = 5 , and /i = 0.1, 0.06, and 0.02 (from top to 
bottom) 
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Fig. 3b 




